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SECTION  ONE 


INTRODUCTION 

The  use  of  fluid  injection  and  combustion  in  the  separated  region  behind 
blunt-based  projectiles  or  vehicles  has  been  a  matter  of  interest  for  some 
time,  ever  since  it  became  known  that  the  base  pressure  could  be  raised  by 
this  means.  Some  of  the  early  work  is  listed  in  References  1  and  2.  However, 
as  results  began  to  be  accumulated,  two  things  became  apparent.  First,  while 
a  small  increase  in  base  pressure  was  relatively  easy  to  obtain  (and  this  ease 
implied  high  specific  impulse),  base  pressures  significantly  above  ambient 
values  were  much  harder,  and  perhaps  impossible,  to  reach.  Second,  at  higher 
flight  Mach  numbers  the  base  drag  is  a  decreasing  fraction  of  the  total  drag. 

Both  these  factors  indicated  a  need  for  higher  base  pressures,  and  methods 
of  overcoming  the  ambient  pressure  limit  began  to  be  considered.  That  search 
is  still  in  progress.  The  main  method  of  approach  is  to  try  different  means 
of  introducing  the  injectant,  primarily  by  means  of  lateral  injection.  Changes 
in  vehicle  and  injector  geometry  and  in  the  injectant  used  have  also  been 
made.  However,  substantial  success  has  not  yet  been  achieved,  as  can  be  seen 
by  considering  the  two  most  extensive  and  most  recent  experimental  programs  in 
this  area. 

3 

In  the  first  of  these  ,  a  strut-mounted  cone-cylinder  was  used  with  lateral 
injection  of  a  liquid  fuel  (triethylaluminum,  or  TEA)  upstream  of  the  base,  at 
Mach  numbers  up  to  seven.  The  experimental  results  are  classified,  but  may  be 
summarized  by  saying  that  base  pressures  (Fg/P^)  significantly  greater  than 
one  were  only  obtained  under  transient  conditions.  The  authors  of  Reference  3 


conjecture  that  with  proper  design  the  transient  results  may  he  sustained, 
hut  this  has  not  yet  been  shown  to  be  true. 
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The  second  program  of  interest  used  as  injectant  the  exhaust  of  a 
number  of  solid  propellants,  most  of  winch  were  substantially  fuel-rich.  Here 
base  pressure  ratios  as  high  as  2  were  obtained,  but  some  questions  of  tunnel 
influence  exist  for  all  base  pressure  ratios  greater  than  one. 

These  results  seem  to  support  the  conclusion  that  substantial  base  pressure 
ratios  (say  P^/P^  -  2)  are  difficult,  but  perhaps  not  impossible,  to  obtain. 

It  is  clear  that  the  potential  range  of  vehicle  and  injection  geometries,  and 
of  injectants,  is  very  wide,  and  only  isolated  points  have  been  tested  to  date. 

We  are  encouraged  to  hope  for  higher  base  pressure  ratios  as  a  result  of 
the  concept  developed  by  Strahle^ .  Here  the  injectant  penetrates  the  surroundinq 
air  stream  and  bums  relatively  far  from  the  body,  in  an  annulus  surrounding 
tiie  wake.  Under  these  conditions  the  flow  passing  within  trie  combustion  annulus 
is  deflected  towards  the  centerline  by  a  shock  about  the  annulus.  In  principle 
rather  high  pressures  can  be  generated  by  such  a  shock,  especially  at  high  Mach 
numbers,  and  such  high  pressures  are  readily  transmitted  to  the  base  throunh 
the  low  velocity  separated  region,  ilo  experimental  confirmation  of  the 
feasibility  of  this  method  is  yet  available.  This  is  probably  due  at  least  in 
part  to  the  fact  that  such  an  annular  heat  release  region  is  difficult,  perhaps 
impossiule,  to  obtain. 

At  this  point,  then,  it  appears  that  experimental  data  are  required  to 
answer  the  question  of  what  performance  can  be  achieved  using  external  burning 
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propulsion.  But  it  seems  likely,  as  illustrated  By  the  existing  data  base, 
that  high  performance  is  available  only  in  a  very  limited  region  of  that 
multi-dimensional  space  defined  by  all  the  many  variables.  If  this  is  true, 
an  experimental  search  for  this  limited  region  can  be  prohibitively  expensive, 
and  an  increase  in  the  performance  of  external  burning  propulsion  (a  term 
used  here  to  denote  any  propulsion  or  drag  reduction  scheme  using  the 
introduction  of  reactive  fluids  in  or  near  the  base)  is  best  approached  using 
an  anlaytical  model,  by  which  the  performance  of  many  different  candidate 
systems  can  oe  calculated.  Such  a  model  would  allow  experiments  to  oe 
directed  towards  system  configurations  of  high  performance  potential. 


SECTION  TWO 


ANALYTICAL  MODEL 


Pre-Existing  Models 

A  model  of  base  burning  propulsion  can  be  developed  by  modifying  existing 

base  flow  methods  to  include  the  effects  of  injection.  Tiie  Crocco-Lees  and 
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Korst-Chapman  approaches  have  been  the  subjects  of  such  modification,  ’ 
but  results  have  not  been  entirely  satisfactory,  due  partly  to  the  non¬ 
trivial  problems  involved  in  treating  axial  symmetry,  and  partly  to  tne 
complexity  of  these  methods,  when  regions  of  different  fluids  must  be 
distinguished. 

Another  approach  is  to  develop  an  analysis  directly  for  the  base  burning 
case.  Tin's  was  done  by  Schetz  et  al .  for  base  bleed  only*7  and  then  extended 

Q 

to  include  lateral  injection.  The  model  used  accounts  for  many  important 
features  of  the  flow,  including  downstream  mixing  and  coupling  of  the  viscous 
wake  flow  with  the  inviscid  surrounding  flow.  The  model  provides  an 
interesting  Insight  into  the  case  of  a  high  Mach  number  flow  of  low  static 

temperature,  but  with  hinh  staonation  teroeratnrp  flow  there  is  a  ouestinn 
whether  chemical  reaction  heats  or  cools  the  flow. 

However,  the  model  of  Schetz  et  al .  has  limitations  that  make  it  less 
useful  for  lateral  injection.  As  it  stands,  it  does  not  treat  the  injection 
process  in  any  detail.  Further,  in  the  near  wake  region,  it  appears  to  treat 


the  cavity  and  the  surrounding  shear  layer  as  a  single  one-dimensional 
(radially  averaged)  flow.  Although  this  is  a  good  assumption  for  massive 
base  bleed,  it  is  probably  not  as  good  for  lateral  injection. 

The  present  work  is  an  attempt  to  improve  the  treatment  of  lateral 
injection  by  usirg  existing,  detailed  computer  models  of  reactive  liquid  and 
gas  injection  (for  the  liquid  injection  case,  see  Ref.  9),  combined  into  a 
code  we  refer  to  as  BBLIP  (for  Base  Burning/Lateral  Injection  Propulsion), 
the  original  version  of  which  was  produced  by  MDAC  under  contract  to  the 
Ballistic  Missile  Defense  Advanced  Technology  Center  (BMDATC)J® 

BBLIP  calculates  the  flow  field  produced  by  a  single  transverse  jet 
upstream  of  the  base,  and  then  at  the  base  plane  combines  these  individual 
jet  flows  into  a  single  annulus  of  injectant,  of  which  the  inner  boundary  is 
the  wake  cavity,  while  the  outer  boundary  is  matched  in  pressure  and  direction 
to  the  surrounding  flow.  The  model  is  entirely  inviscid,  but  the  presence 
of  the  injectant  annulus  avoids  the  indeterminacy  which  has  long  been  known 
to  exist  in  a  single-fluid  inviscid  analysis;  in  fact,  the  problem  is  mathe¬ 
matically  closed  without  requiring  consideration  of  the  reattachment  pressure 
rise,  and  so  this  can  in  principle  be  calculated. 

The  individual-jet  analysis  was  originally  developed  for  fuel-rich  jet 

11-13 

interaction  control  systems  .  It  considers  two  adjacent  gas  streams 
calculated  stepwise  downstream,  coupled  by  the  requirement  of  pressure  equality 
across  the  dividing  steamline.  The  inner  stream  includes  the  jet  flow  plus 
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all  the  air  that  reacts  with  it;  the  outer  stream  is  air  only.  The  region  neai 
the  orifice,  where  both  streams  contain  large  velocity  components  normal  to 
the  wall,  is  given  special  consideration.  Mixing  is  treated  semi-empirically , 
using  the  limited  data  on  concentration  profiles  in  supersonic  flow  appearing 
in  the  literature.  Equilibrium  chemistry  is  used  to  treat  reaction,  and  any 
proportion  of  the  inner  flow  may  be  assumed  to  react.  An  approximate  method 
is  used  for  calculating  the  entire  flow  about  the  jet,  both  for  convenience 
and  because  the  accuracy  of  the  method  is  appropriate  to  other  parts  of  the 
analysis. 

The  flow  field  as  conceptually  divided  up  for  analysis  is  siiown  in 
Figure  1.  The  analysis  proceeds  as  follows.  Calculation  of  tiie  local  flow 
(region  I)  yields  the  state  variables  describing  the  undisturbed  air  stream 
at  tne  point  of  injection.  These  state  variables  are  then  assumed  to  ue 
uniform  everywhere  upstream  of  tne  Shock  around  the  jet.  Tney  are  used,  . it 
tne  state  variables  of  tne  jet,  to  calculate  tne  size  of  t.he  obstacle,  assumed 
shaped  like  a  quarter  sphere,  equivalent  to  that  presented  by  the  jet  lo 
the  external  flow.  This  calculation  allows  an  entropy  increase  due  to  tne 
complex  shock  system  present  in  the  jet,  and  includes  the  effects  of  jet 
and  free  stream  liach  numbers  and  jet  cant  angle, 

Tne  flow  of  tne  external  stream  around  this  obstacle  is  obtained  from 
the  analysis  of  region  II,  which  assumes  inviscid  compressible  flow  around 
a  quarter  sphere.  This  calculation  yields  both  tne  siiock  snape  and  tne 
pressure  along  tne  dividing  streamline,  the  latter  needed  for  tne  momentum 
oa lance  in  region  III,  to  follow. 


In  region  III,  a  mixing  analysis  is  used  to  define  a  point  above  wiiicn 
tiie  fraction  of  jet  fluid  is  small.  Tin's  cutoff  value,  different  values  of 
which  can  be  used  in  regions  III  and  IV,  is  experimentally  determined.  Toon, 
based  on  the  flow  field  computed  in  region  II,  a  streamline  (ttie  dividing 
streamline)  is  followed  back  from  this  point  to  the  shock  and  into  the 
undisturbed  flow  upstream. 

Region  III  is  calculated  one-dirnensional ly,  invoking  conservation  of 
mass  and  momentum.  Radial  gradients  and  vorticity  due  to  turning  the  jet, 
both  known  to  occur,  are  neglected.  On  the  basis  of  experimental  results, 
tiie  flow  in  region  III  is  assumed  frozen. 

Region  IV  is  uuilt  up  of  increments,  in  each  of  whicn  the  calculation 
is  the  same  and  similar  to  that  of  region  III:  a  one-dimensional  analysis 
eased  on  the  conservation  equations  and  accounting  for  heat  release  by 
equilibrium  chemistry,  matching  pressure  with  the  outer  flow  (region  V) 
at  the  dividing  streamline,  and  neglecting  viscosity,  vorticity,  and  normal 
gradients . 

In  region  V  the  outer  flow  and  shock  are  calculated  stepwise  downstream, 
at  each  step  iterating  with  region  IV  on  the  slope  of  tne  dividing  streamline, 
requiring  pressures  to  match  across  it.  The  outer  flow  is  assumed  inviscid, 
compressible,  and  axially  symmetric  about  a  streamwise  line  in  tne  vehicle 
surface. 


At  the  base  plane,  the  individual  nozzle  flows  are  transformed  to  a  single 
annulus  concentric  with  the  body  and  axially  symmetric,  and  the  inner  boundary 
of  this  annulus  turns  inward  as  shown  in  Figure  2  to  an  angle  In  this 
early  model,  the  inner  boundary  pressure  is  calculated  using  a  simplified 
recirculation  model,  which  includes  a  term  accounting  for  the  addition  of  base 
bleed,  and  thus  allows  evaluation  of  the  possible  synergistic  effect  of  com¬ 
bining  base  bleed  and  lateral  injection.  Tiie  outer  boundary  pressure  is 
calculated  from  expansion  of  the  outer  flow  through  an  angle  0c  +  ow  -  6, 
where  is  the  cone  half  angle  and  6  is  the  divergence  angle  between  inner  and 
outer  annulus  boundaries,  dote  tnat  because  of  axial  symmetry,  in  general 
d  >  0  even  if  the  annulus  flow  is  taken  as  frozen  and  without  mixing.  The 
streamwise  development  of  the  injectant  annulus  is  calculated  using 
equilibrium  chemistry,  in  zones  of  varying  concentration,  to  prevent  losing, 
via  cross-stream  averaging,  the  highly  nonlinear  effect  of  concentration  on 
heat  release.  Conservation  of  mass  and  momentum  in  the  injectant  annulus 
also  is  required. 

Entrainment  of  air  is  accounted  for  via  the  angle  a  (Figure  2),  which  is 
considered  to  be  unknown,  to  be  determined  by  comparing  program  predictions 
with  the  results  of  experiments.  The  external  shear  layer  somewhat  resemules 

the  shear  layer  bounding  a  coaxial  jet,  for  wiiicli  data  exist  that  could  be 
used  to  calculate  a.  It  is  likely,  however,  that  tiie  mixing  half-angle  of 
the  external  snear  layer  will  depend  on  the  geometry  at  the  injectors  wnere 
that  shear  layer  originates,  and  therefore  o  is  preserved  as  an  unknown. 
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OUTER  BOUNDARY 


ure  2.  Injectant  Annulus  for  Base  Burning  Propulsion  Using  Lateral  Injection;  Geometry 
Approximately  to  Scale 


In  calculating  base  pressure,  normal  pressure  gradients  in  the  injectant 
annulus  are  assumed  negligible.  The  calculation  advances  as  follows:  First, 
bw  is  chosen.  Second,  <s  is  varied  until  Pexit  "  ^  ^Pouter  "  Pinner^ 

Pouter  meets  a  Prese*  criterion.  The  resulting  pressure  is  the  desired 
solution. 

The  calculations  described  above  are  also  outlined  in  Reference  10,  which 
contains  some  references  to  other  work.  A  more  complete  description,  together 
with  an  experimental  program  and  the  use  of  the  results  to  determine  a, 
can  be  found  in  Reference  14. 

Modifications  to  the  BBLIP 

At  the  start  of  the  present  work,  BBLIP,  as  described  above  and  in 

References  10  and  14,  was  deficient  in  several  respects.  First,  it  was 

developed  to  run  at  high  Mach  numbers  only.  This  sufficed  for  the  interceptor 

application,  but  the  present  context  required  the  ability  to  run  below  Mach  5, 

down  to  Mach  2  and  perhaps  below.  Second,  further  consideration  of  the 

reattachment  process,  in  light  of  the  r”^  terms  in  the  axi symmetric 

equations  of  motion  led  to  a  reformulation  of  the  model  and  a  return  to  the 
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physics  observed  by  Chapman.  And  finally,  improvements  were  needed  in 
accounting  for  axial  symmetry  in  both  the  injectant  annulus  and  the  external 
flow.  These  changes  are  described  in  detail  in  the  rest  of  this  section. 

Low  Mach  Humber 

At  the  start  of  the  present  program,  individual-jet  calculations  upstream 
of  the  base  could  not  be  made  below  about  M^  =  b.  Two  changes  were  required 
in  the  individual-jet  calculation,  as  documented  in  Reference  11-13,  in  order 
to  operate  in  tne  range  2  -  M^  -  5. 


Upstream  region:  The  original  calculation  in  the  upstream  region 

(defined  in  Figure  1)  treated  inviscid  flow  around  a  quarter  sphere  to  yield 

two  required  results:  (1)  the  jet  shock  shape,  and  (2)  pressure  along 

the  dividing  streamline  in  the  upstream  region.  The  latter  is  used  in  the 

PdA  term  in  the  inner  flow  analysis  (region  III  of  Figure  1).  The  radius 

of  the  quarter  sphere  was  a  scale  length  called  (not  very  accurately)  the 

jet  penetration  height  h,  calculated  from  local  stream  and  jet  conditions 

by  a  method  due  to  Kail  is. ^  The  flow  was  calculated  by  the  approximate 
17  18 

method  of  Maslen  as  extended  by  Jackson;  the  resulting  shock  shapes  were 
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found  to  compare  reasonably  well  with  high-  ML  experiments.  * 

Adjustments  to  the  computer  program  made  calculations  possible  at  Mach 

numbers  as  low  as  2.42,  but  this  did  not  meet  the  present  requirements.  A 

21 

method  was  therefore  developed  which  uses  empirical  shock  shapes  but  retains 
17  18 

the  Maslen  -  Jackson  integrations  inward  from  the  shock  through  the 
shock  layer  to  the  dividing  streamline. 

Results  of  this  method  were  compared  with  results  of  the  original 
Maslen- Jackson  method,  as  embodied  in  a  code  known  as  UASJAC.  The  latter  had 
been  shown  to  compare  reasonably  with  experiment.  In  the  range  -  2.42  where 
both  are  operational. 

In  making  these  comparisons,  it  became  clear  that  the  values  of  pressure, 
velocity  and  entropy  resulting  from  the  empirical-shock  calculation  differed 
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substantially  from  those  calculated  uy  IIASJAl,  under  conditions  where  notn 
calculations  can  be  made.  Since  MASJAC  nas  been  compared  with  experiment, 
and  since  ttie  empirical  expressions  were  developed  comparing  shock  shapes 
ratner  than  internal  properties  (tne  latter  oeinq  much  more  sensitive  than 
the  former),  it  was  decided  to  adjust  ttie  empirical  expressions  to  natch  the 
ilASJAC-calculated  properties.  Two  parameters  in  the  empirical  expressions 
can  be  simply  adjusted:  standoff  distance  and  nose  radius  of  curvature.  Tne 
importance  of  ttie  former  is  trivial  in  the  present  context,  and  so  only  channes 
in  the  nose  radius  of  curvature  were  considered.  A  value  of  II,.  =  'J.b  x 
(empirical  value)  gives  good  agreement  with  iiASJAC  over  the  range 
<?.4  <  fi^  <  9.8.  This  expression  is  therefore  used  in  utJLIP  for  all  values  of 

‘i- 

Downstream  Region:  At  each  downstream  increment  in  ttie  individual-jet 
calculation  (see  Figure  1)  the  jet  shock  and  inner  flow  are  calculated.  Tne 
shocklayer  calculation  again  is  done  using  1-iASJAC.  At  lower  values  of  liL,  the 
jet  stiock  moves  away  from  the  inner  flow,  and  the  computational  rays,  wnicti 
are  shown  in  Figure  1  as  lines  normal  to  ttie  jet  shock,  become  longer.  At 
some  point,  the  possibility  exists  that  these  lines  may  cross  as  the  shock 
angle  is  adusted  during  the  process  of  reaching  a  solution. 

When  this  happens,  problems  can  usually  be  avoided  by  choosing  a  smaller 
number  of  increments.  A  range  of  dov/nstream  increment  number  ( L>i  1 1 )  has  ween 
found  allowing  runs  to  be  made  at  all  conditions.  Tnese  results  are  sriown 


in  figure  3,  which  shows  that  for  any  downstream  region  more  than  1/2  and 
less  tiian  about  lu  scale  heights  long,  a  range  of  Dill  exists  which  allows 
the  program  to  run.  Values  too  high  or  too  low  incur  problems  which  are 
not  theoretical  but  computational. 


Since  UNI  must  be  an  integer,  the  following  expression,  rouqnly 
corresponding  to  the  line  in  Figure  3,  is  used  in  UbLIP: 

>  6zds 

UNI  =  Integer  -  1.4  —r — 

one  expects  UNI  to  affect  tne  answer,  in  tnat  mesh  size  generally  has  sucn 
an  effect  on  numerical  analysis.  Tiiis  effect  is  shown  in  Figure  4,  in  wi.icn 
the  usual  trend  of  reaching  an  asymptote  with  decreasing  mesh  size  (here 
increasing  UNI)  is  shown,  together  witu  the  negligible  magnitude  of  tnis 
change . 


Wake  Closure 

It  tias  been  known  since  Chapman's  work J J  in  19ol,  tnat  a  conical  base 
cavity  cannot  be  truly  at  constant  pressure.  Chapman  also  observed,  from 
shadowgraphs  of  the  near  wake,  that  tne  wake  looks  like  a  cone  followed  > 
a  cylinder  of  radius  R^.  Examples  of  this  geometry  can  be  seen  in  Chapman's 
report* J  and  on  page  349  of  Liepmann  and  Roshko.^  This  shape  is  formed,  not 
by  streamlines,  out  oy  the  outer  boundary  of  the  turbulent  region,  however, 
in  supersonic  flow  the  mixing  half-angle  is  small  (as  will  be  discussed  below), 
and  tnus  tiie  outer  boundary  of  the  turbulent  region  approximates  a  streamline. 
It  was  therefore  decided  to  change  the  wake  model  (which  up  to  now  carried 
tne  inviscid  streamlines  down  to  Rw  =  U,  as  shown  in  Figure  2)  to  require 
it  to  end  at  non-zero  R^. 
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(NOTE:  OTHER  FLOW 
PROPERTIES  ARE  LESS 
SENSITIVE  TO  DNI 
THAN  VELOCITY: 


Figure  b  shows  Chapman's  data,  together  with  a  correlation  developed  here. 
It  is  not  entirely  clear  whether  to  use  the  injectant  annulus  Mach  number 
or  the  local  tlach  number  fi^,  in  the  correlation.  Arrangements  for  either 
can  be  made.  At  present  is  used. 

Structure  of  the  Present  Model 

The  present  model  uses  the  individual-jet  model  to  calculate  injectant 
conditions  at  the  base  plane,  and  all  these  jet  flows  are  mapped  into  an 
annulus  there,  just  as  in  the  previous  model.10  The  downstream  evolution  of 
tiiis  injectant  annulus  is  now  calculated.  Instead  of  making  this  calculation 
in  a  number  of  steps,  as  is  done  upstream  of  the  base  plane,  a  single  step 
is  taken  from  the  base  plane  to  the  point  where  the  annulus  inner  radius 
equals  Ky .  This  single  step  uses  information  about  the  annulus  geometry,  as 
shown  in  Figure  6.  The  calculation  uses  conservation  of  mass  and  of  axial 
momentum,  and  thermochemistry. 

Conservation  Equations 

The  conservation  equations  become 

piA1  *  pi„  "  <rB2  -  R»2>  +  Pout  u  >  <RtUrn2  *  <RB  *  Rl>2> 

*  Pout  d  ’  (Rturn2  '  R32>  +  ”l  vi  +  “u  V1  +  S"'d  v2  (1> 

-  (m^  +  6mu  +  6md)  Vi+^  =  P,+1  n  (R^  -  Rw  ) 
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Figure  6.  Injectant  Annulus  Geometry  (Approximately  to  Scale  for  =  3,  V^Ajj  =  0.0313 


Finding  the  latter  8  parameters  is  the  key  to  the  problem.  In  FORTRAN 
these  are  respectively  (POUTU,  POUTD,  PIN,  RTURN,  DMDTU ,  DMDTD ,  VI,  V2). 


f  i_nd ing  the  Required  Parameters 

RTURN:  The  injectant  annulus  expands  at  the  vehicle  base  plane  through  the 
angle  (0c  =  0W).  In  reality  this  expansion  takes  place  in  a  fan-shaped  region 
centered  at  point  A  in  Figure  6.  For  convenience  we  assume  that  it  can  be 
represented  as  occurring  at  a  single  ray,  OA  in  Figure  6,  placed  at  an  angle 
,i  to  the  upstream  flow.  The  value  of  a  is  found  by  requiring  the  Mach  number 
change  through  the  expansion  to  satisfy  continuity.  Then,  since  the  expansion 
is  assumed  isentropic,  all  other  values  also  change  correctly,  and 


v  ,  nii  +  “» *  «n _ 

1+1  n  <R32  -  rm2) 


Substituting  (2)  in  (1)  yields  a  quadratic  for  R^  ,  which  requires  values 
of  the  following  parameters: 

( P ,  A^,  rru ,  v.,  Rg,  R-j)  obtained  from  initial  conditions 


(Rw)  obtained  from  experiment 


(RT n-+i )  obtained  from  thermochemistry 


(P.j+l)  obtained  from  assumption  P.+^  =  PQut  ^  =  P.^ 


(Pout  u’  Pout  d*  Pin ’  Rturn’  fV  «"e’  V  V  to  be  stained. 


— 1 


The  geometry  is  shown  in  Figure  7.  Using  continuity  and  P  =  p RT ,  we  have 


P,  R, 


M1  yfh 


-  ?2  r2  ^2 


i1 


which  yields 


1  +  ^  »i2 

1  +  ^  m22 


A1  so 


rl  r2 

sin  a  =  ,  sin  (a  +  0)  =  j-=- 


(3) 


(4) 


where  L  is  the  length  of  the  equivalent  single  ray  (Figure  7).  Solving 
(3)  and  (4)  for  a  yields 


tan 


-1 


sin  0 


1  *  *r 


Y  +  1 

27^7 


-  cos  0 


(5) 


Since  we  are  discussing  flow  in  the  annulus,  M-j  =  Mi  and  is  obtained  from 
a  Prandtl-Meyer  expansion  through  0;  Figure  6  shows  0  =  0w  +  Qc.  Thus 
interpreted,  (5)  yields  a  from  known  initial  conditions  and  an  assumed  value 


From  Figure  6,  a  determines  Rturn-  Assuming  small  angles, 

0  +  <f>  +  ex 

Rtll„  =  Rn  +  R,  — - 

turn  B  1  a  -  a 


(6) 
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Thus  the  first  of  the  required  variables  is  determined,  based  on  knowledge 


of  the  angles  in  Figure  6,  particularly  G  . 

W 

Pressures:  Assuming  the  angles  known,  the  bounding  pressures  are  easily 
determined  if  Figure  6  describes  planar  flow.  In  that  case  pressures  are 
constant  along  straight  streamlines,  ,>Ogtu  =  ’'i-  P1„  ■  P1  *  SPPH  (0C  * 

and  pout  d  '  ’’out  u  *  sppn(0c  *  0W  '  wBere  the  notation  6PpH  (6) 
refers  to  the  pressure  drop  in  a  Prandtl-Meyer  expansion  through  the  angle  G. 


However,  Figure  6  represents  not  planar  but  axisymmetric  flow.  In  such  a 
flow  the  pressure  changes  along  any  straight  streamline  that  is  not  parallel  to 
the  axis.  This  is  readily  calculated  using  the  method  of  characteristics, 
but  in  the  present  analysis  the  flow  field  of  Figure  6  must  be  calculated 
many  times,  and  an  approximate  method  is  needed. 


,  23  24 

Such  a  method  was  developed  by  Webb  and  used  by  Mehta  and  Strahle  . 

x'  This  is  based  on  Webb's  observation  that  the  effects  of  streamline  convergence, 
in  a  typical  near-wake  geometry,  are  communicated  mainly  by  the  right-running 
characteristics,  the  "reflected"  left-running  characteristics  playing  a 
relatively  minor  role  (see  Figure  8).  This  approximation  then  allows  us  to 
use  the  characteristic  relation  along  the  right-running  characteristic,  which 
in  finite  difference  form  is 


sin  0  sin  ug  r3  “  r2 
sin  (0  +  ^2 


(7) 


to  give  the  change  in  Prandtl -Meyer  angle  v  as  a  function  of  the  changes 
in  flow  direction  0  and  radius  r,  and  the  Mach  angle  y.  The  numbering  and 
location  of  points  is  shown  in  Figure  8. 


Note  that  (7),  applied  to  the  cavity,  allows  the  use  of  a  constant  pressure 
boundary  but  then  results  in  non-constant  direction  (the  same  result  Chapman15 
obtained).  However,  on  the  basis  of  observed  geometry16,  we  take  the  other 
option,  assuming  direction  constant  (i.e.,  e3  -  ©2  =  0)  and  finding  v(r) 
and  thus  P(r). 


The  term  (r3  -  in  (7)  is  obtained  as  follows:  from  Figure  8,  it 

is  clear  the  following  expressions  hold: 

L03  =  z3/cos  9 

z3  =  (ro  -  r3)  ctn  0 

1-23  =  1/2  LQ3/cos  p 

rZ  =  r3  +  1-23  cos  ^2  ~  9  “  u) 


From  these  we  obtain 


r3  ~  r2  _  R  -  1 

r2  1  +  R  (|  -  1 ) 

where 


A  = 


sin  (9  +  p) 
sin  9  cos  p 


(8) 


(9) 
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and  R  =  radius,  nondimensionalized  in  such  a  way  that  the  final  value  is  1; 


thus  for  the  inner  flow  R  =  r^/RW;  for  the  outer  flow  upstream  of  R^urn> 

R  =  r3^turn’  ^or  tlie  downstream  outer  ^low,  R  =  r-j/R-j. 

Thus  (7),  with  (8)  and  (9),  gives  v(r),  and  the  Prandtl-Meyer  relation 
gives  P(r).  Then  average  pressures  for  use  in  the  momentum  equation  are 
obtained  as  follows: 


1 


(10) 


POUTU  = 


turn 


-  K +  R 


•)2li 


f 

J  Rn  + 


2tt  r  dr  P(r) 


(ID 


turn 


POUTD  = 


_ 1 

'  (R32  '  RtUrn2) 


r  dr  P(r) 


(12) 


Entrainment:  The  injectant  annulus  entrains  air  from  the  external  flow  via 
mixing  in  the  outer  shear  layer.  This  mixing  is  accounted  for  by  inclining 
the  outer  boundary  of  the  annulus  at  an  angle  a  to  the  local  streamlines;  o 
is  the  mixing  half-angle  defined  by  the  maximum  slope  of  the  concentration 
profile. 


I 
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(Its  relation  to  the  similarity  length  scale,  also  often  denoted  o,  is 
discussed  by  Brown  and  Roshko.  )  The  area  through  which  the  entrained  flow 
passes  is  thus  defined,  but  before  its  mass  flow  rate  can  be  obtained,  its 
density  and  velocity  must  be  known. 

These  are  found  from  the  average  pressures,  by  assuming  the  flow  isen- 
tropic.  The  following  relations  are  used: 

Y  - 
2 


In  the  above,  and  are  values  outside  the  annulus  at  the  base  plane,  wnile 
(and  therefore  M2)  are  calculated  from  P  (Rturn)»  obtained  the  same  way  as 
P(r)  in  (10)  -  (12),  and  then  expanded  around  the  comer  at  Rturn. 

Mixing  Half-Angle  o:  Ideally,  o  would  be  determined  by  analysis  of  the 
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external  shear  layer.  In  1972,  when  the  Langley  Conference  on  free  turbulent 
shear  flows  was  held,  it  was  known  from  experiment  that  for  a  turbulent  shear 
layer,  above  a  stagnant  region,  mixing  is  decreased  from  the  subsonic  case 

14 

when  the  flow  is  supersonic.  In  fact,  a  reasonably  good  fit  to  the  data  is 

o  -  - - -  (17) 

1  +  M2 

The  ability  to  calculate  this  sort  of  behavior,  however,  is  another  matter. 

Unly  a  few  analyses  show  anything  like  this  Mach  number  dependence,  and  though 
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more  work  has  been  done  on  the  subject  since  then,  a  reliable  prediction  of 
o  does  not  yet  exist  even  for  this  relatively  simple  case.  The  present  appli¬ 
cation  is  further  complicated  by  including  supersonic  flows  of  different 
gases  on  either  side  of  the  shear  layer,  vortical  structures  from  tne  indi¬ 
vidual  jets,  and  streanwise  and  transverse  pressure  gradients.  Here  we  evade 
these  difficulties  by  leaving  o  to  be  determined. 


28 


SECTION  THREE 


NUMERICAL  ASPECTS 

Method  of  Solution 


The  unknowns  in  (1)  and  (2)  have  now  been  obtained,  leaving  one  equation 
for  R3,  the  annulus  outer  radius  at  the  cone-cylinder  transition,  but 
requiring  the  angles  shown  in  Figure  6.  The  method  of  solution  is  as  follows: 

(1)  A  value  of  6,  the  annulus  divergence  angle,  is  assumed. 

(2)  A  value  of  0W  is  assumed. 

(3)  From  these  and  the  base  plane  values,  PIN  and  POUTO  are  calculated. 

(4)  Adjust  0  ,  at  constant  6,  until  c  =  (PIN-P0UTu)/P0UTD  is 

w  p 

sufficiently  small. 

(t>)  Set  P^+i  =  (PIN  +  P0UTU)/2  and  calculate  R3  from  equations  (1)  and 

(2). 

(6)  Calculate  R3S  from  the  geometry  shown  in  Figure  6. 

(7)  Adjust  6  until  cR.j  3  (R3  -  R3S)/R3  is  sufficiently  small. 

A  flow  chart  of  the  resulting  computer  program  is  shown  in  Figure  9. 

Solution  Behavior 

Multiple  Solutions 

The  formulation  required  for  even  the  highly  approximate  version  of  reality 
described  in  the  previous  paragraphs  is,  as  can  be  seen  above,  conceptually 
simple  but  mechanically  fairly  complex.  Difficulties  have  been  encountered 
in  obtaining  solutions  under  some  sets  of  conditions.  In  this  section  we 
discuss  some  of  these  problems. 


Whatever  the  cause,  the  problem  most  often  seen  shows  up  as  multiple 
solutions,  or  as  a  failure  to  reach  a  solution.  Referring  to  the  simplified 
solution  outline  presented  above,  wfiat  often  happens  is  the  following.  Steps 
(1)  -  (4)  can  usually  be  completed  without  difficulty,  as  can  steps  (S)  and 
(6),  which  are  purely  mechanical.  But  either  no  6,  or  more  than  one,  can 
sometimes  be  found  to  satisfy  step  (7). 

how  it  seems  clear  that  neither  no  solution,  nor  prouably  multiple 

solutions,  can  correspond  to  physical  reality.  However,  this  difficulty  is 

believed  to  be  related  to  the  physics,  possibly  as  follows.  Consider  PIN  and 

POUTd,  which  are  forced  to  match  by  adjusting  0  and  <$.  (It  is  found  that 

0  and  6  are  linearly  related,  so  it  suffices  in  this  discussion  to  talk 
w 

in  terms  of  one  of  them,  say  0  .)  In  this  way  the  geometry  is  obtained,  which 

V* 

then  yields  an  exit  radius  R3S  and  allows  calculation  from  the  conservation 
equations  of  another,  R3,  which  must  match  it. 

But  PIN  and  POUTD  are  decreasing  functions  of  Gw.  Both  result  from 
expanding  different  flows  (respectively  the  injectant  and  external  flows), 
and  often  from  nearly  the  same  value.  (These  initial  values  will  be  the 
same  if  the  annulus  streamlines  upstream  of  the  turning  ray  or  expansion  are 
parallel  to  the  centerline,  because  by  the  upstream  individual-jet 
calculation,  inner  and  outer-flow  pressures  are  forced  to  be  equal  at  the  base 
plane.)  if  the  expansion  begins  at  the  same  value,  and  if  the  functions  P(e^) 
are  similar,  no  solution  can  exist  (Fig.  Id).  If  the  functions  are  dissimilar, 
one,  two,  or  no  solutions  can  exist,  as  that  Figure  snows. 
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Figure  10.  Examples  of  Possible  Solution  Behavior 


Different  behavior  would  be  found  if  one  of  the  two  functions  were  an 
increasing  function  of  8  .  As  it  is,  we  may  expect  difficulties  with  stability 
as  well  as  number  of  solutions;  Figure  10(c)  shows  an  example  where  the 
solution(s)  would  change  substantially  with  a  very  small  change  in  initial 
conditions  or  functional  dependence. 

An  example  of  multiplicity  or  absence  of  solutions  is  shown  in  Figure  11, 
which  shows  eR3  as  a  function  of  <S.  As  shown,  single  solutions  seem  to  exist 
for  M  =■  2  and  M  *  5,  and  no  solutions  for  M  =  2.5,  3  and  3.5.  However, 

Figure  11  (as  noted  there)  is  for  a  =  0.6  deg.,  and  it  has  been  found  that 
varying  a  has  the  effect  of  translating  the  whole  pattern  up  or  down  (increasing 
o  moves  the  pattern  oown)  without  changing  its  character.  In  particular,  an 
increase  in  a  might  bring  the  pattern  far  enough  down  that  line  b  in  Figure  11 
would  coincide  with  eR3  =0.  If  that  were  the  case.  Figure  11  would  show  a 
single  solution  for  M  *  2,  multiple  solutions  for  M  =  3  and  3.5,  and  probable 
multiple  solutions  for  M  =  2.5  and  5. 

These  solutions  can  be  denoted  plus  and  minus  solutions,  based  on  the 
sign  of  deR3/d6.  Thus  the  single  solution  for  M  =  2  is  a  plus  solution,  and 
that  for  M  =  5  is  a  minus  solution.  In  fact  the  above  discussion  implies  that 
we  can  often,  perhaps  always,  adjust  a  to  obtain  two  solution  branches  over 
much  of  the  Mach  number  range.  The  "correct"  branch  might  then  be  chosen  by 
some  as  yet  undetermined  means. 

But  when  this  is  done,  two  observations  are  made.  First,  only  the  plus 
solution  can  be  reached  at  Mach  2,  and  only  the  minus  at  Mach  5.  A  jump  in 
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base  pressure  results  on  transition  from  one  solution  to  the  other  and  as 
ooth  solutions  exist  over  a  wide  range,  transition  is  not  uniquely  located. 
Second,  neither  solution  really  behaved  right,  both  in  terms  of  P(M)  and 
of  0W  (M) .  Observation  of  shadowgraphs  for  a  wide  range  of  conditions 
showed  10  <  6W  <  15  degrees.  This  is  without  injection,  and  yet  it  is 
believed  that  the  same  range  should  be  seen  with  injection,  at  least  in  the 
absence  of  very  large  heat  release.  Neither  solution  branch  stays  in  tne 
correct  6W  range  over  the  whole  Mach  number  range. 

Another  option  exists,  typified  by  the  dashed  line  lauded  "C"  in 
Figure  11.  A  single  solution  can  be  made  to  exist  over  the  whole  range  of 
Mach  numbers,  by  adjusting  the  pattern  shown  in  Figure  11  so  that  the  curve 
with  the  highest  minimum  (in  Figure  11  it  is  the  one  corresponding  to  Mach  3) 
to  be  tangent  to  the  “  0  line.  We  now  have  one  solution,  and  only  one, 

at  the  corresponding  Mach  number;  and  this  means  there  is  one  continuous 
solution  over  the  whole  Mach  number  range.  The  tangency  condition  has  in 
fact  produced  a  transition  between  plus  and  minus  solutions  at  the  point  of 
tangency. 

As  will  be  seen  when  results  are  disclosed,  this  approach  is  not  without 
attractions.  It  produces  reasonable  answers  for  both  base  pressure  and  wake 
angle  over  the  whole  Mach  number  range  (although  it  will  be  seen  that  not  all 
the  trends  predicted  are  so  reasonable).  Before  looking  at  results,  however, 
two  questions  need  to  be  discussed. 


« 
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First,  though  this  approach  appears  expedient,  it  is  justifiable?  Is 
the  answer  thus  obtained  the  correct  answer?  The  answer  to  this  question 
is  not  known.  Mathematically  speaking,  no  proof  of  uniqueness  is  advanced 
here.  Considering  the  approximations  involved  and  the  nonlinearity  of 
the  resulting  formulation,  such  a  proof  may  be  impossible  to  obtain  and 
might,  if  obtained,  be  meaningless. 

Physically,  at  each  set  of  vehicle,  flight  and  injection  conditions,  a 

solution  must  exist.  More  than  one  may  exist,  or  transitions  may  exist. 

Although  there  is  some  evidence  suggesting  both  can  occur  under  at  least 
28 

some  conditions,  it  seems  probable  that  base  pressure  or  any  geometrical 
variable  is  a  continuous  function  of  any  of  the  independent  variables. 

Insofar  as  the  approach  under  discussion  produces  continuous  results,  it 
appears  justified  in  this  limited  sense. 

Second,  how  can  the  pattern  be  adjusted?  It  is  known  that  this  can  oe 

done  by  changing  the  mixing  half-angle  a.  However,  o  cannot  be  changed 

freely.  There  is  evidence  that  o  depends  on  the  Mach  numbers  and  densities 

25 

of  the  flows  bounding  the  shear  layer  (see  also  the  discussion  in  Section  II 

11  12  13 

above),  and  that  a  decreases  downstream  of  individual  lateral  jets,  *  * 

presumably  as  the  turning  vortices  die  away.  Other  changes  in  o  cannot  be 
justified,  and  if  they  are  required  in  order  to  produce  continuous  results, 
the  entire  method  comes  under  suspicion. 

It  is  believed  that  there  is  another  way  of  adjusting  the  pattern:  by 
changing  g>b ,  the  external,  upstream  divergence  angle  (Figure  b).  At  present 


this  change  is  only  possible  in  principle,  as  ^  is  included  in  the  external 
flow  formulation  in  such  a  way  as  to  preclude  negative  values,  and  is  not 
included  in  the  annulus  flow  formulation  at  all.  However,  it  is  believed 
that  this  Is  an  area  where  continued  work  might  be  fruitful. 

Radial  Momentum 

Conservation  of  axial  momentum  is  required,  as  described  above  (Equation 
(1)),  but  radial  momentum  is  left  to  fend  for  itself.  The  change  in  pressures 
as  streamlines  change  direction  should  account  for  radial  momentum  conservation. 
This  is  true,  for  example,  in  a  Prandtl-Meyer  expansion  in  planar  flow,  or  in 
an  axisymmetric  flow  calculated  by  the  method  of  characteristics.  In  the 
present  case  an  approximation  to  the  method  of  characteristics  is  used,  and 
the  amount  of  error  involved  in  the  approximation  is  not  known.  The  degree  to 
which  radial  momentum  is  conserved  gives  an  estimate  of  the  errors  involved. 

This  is  done  as  follows.  The  component  terms  of  the  radial  momentum 
equation  are  calculated: 

PAROU  =  POUTU  *  A(r,  z,  9) 

2 

where  A(r,  z,  9)  =  2irrz  +  irz  tan  9 
and  r  =  RBRI 
z  =  LTURN 

e  =  <t>c  +  'J’b  +  ° 

MV ROD  =  -DMDTUP  *  VOUTUP  *  sin  (9C  +  *b) 

PAROD  =  POUTD  *  A(R3S,  HL-LTURN,  9 w  -  6  -  a)  (18) 

MV ROD  =  DMDTDN  *  VOUTDN  *  sin  (9W  -  6) 

MVREX  =  -  mi+1,  vi+1,  sin  (6W  -  ) 
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PAR  IN 


-PIN  *  A(RW,  HL,  0w) 

MVRINL  =  -  m.  sin  (0c  +  <t>b)  (18  cont.) 


Then,  to  estimate  the  degree  by  which  we  fail  to  conserve  angular  momentum, 
we  calculate 


RMOMCK 


algebraic  sum  of  above  terms 
sum  of  absolute  values  of  above  terms 


(19) 


Values  of  RMDMCK  at  solutions  range  from  0.17  to  0.05,  with  no  trend 
evident  from  which  systematic  errors  can  be  deduced.  The  values  seem 
reasonably  small,  and  not  inconsistent  in  absolute  value  with  the  approxi¬ 
mations  made  in  the  analysis  and  described  above.  The  fact  that  only  positive 
values  of  RMOMCK  are  seen  seems  to  indicate  that  the  negative  terms  in  the 
numerator  of  (19)  are  too  small.  Since  PARIN  dominates  the  negative  terms, 
conservation  of  radial  momentum  would  be  well  approximated  by  multiplying 
PIN  —  and  therefore  PB  —  by  a  factor  (1  +  RMOMCK).  This  cannot  be  done  so 
simply,  as  it  would  upset  the  other  conservation  equations  that  have  already 
been  satisfied;  nevertheless  it  is  probably  an  indication  that  the  base 
pressure  we  calculate  is  slightly  too  low. 


SECTION  FOUR 


COMPARISON  WITH  DATA 

Only  one  set  of  unclassified  data  exist  for  gaseous  injectant  and  axisym- 
metric  geometry.4  The  data  labeled  ARC  on  Figure  12  include  only  those  low- 
altitude,  low-M^  points  that  are  claimed  by  the  authors  of  Reference  4  to 
be  free  of  tunnel  interference  effects.  This  conclusion  may  be  questioned 
because  it  appears  that  a  reflected  jet  shock  would  cross  the  wake  axis  upstream 

of  reattachment,  and  because  their  measured  lip  pressures  appear  too  high. 

Some  question  also  exists  Decause  P../P  appears  independent  of  m.  over  a 

tJ  CO  J 

substantial  range.  Nevertheless,  it  is  the  only  unclassified  result 
available. 

Considerable  data  exist  for  base  pressure  without  injection.  These  are 

•  II 

shown  in  Figure  12  as  the  "  m.  =  0  data  range.  These  data  show  some 
evidence  of  an  inflection  point  similar  to  that  shown  between  Mach  J  and 
4  by  the  present  calculations,  and  the  general  shape  is  not  dissimilar  to 
that  resulting  from  the  present  calculations. 

In  Figure  12  these  data  are  compared  with  results  of  the  present  model. 

The  comparison  is  not  exact,  in  the  sense  that  somewhat  different  conditions 
were  used  in  the  calculation  than  were  present  in  the  experiment.  Conditions 
used  in  the  calculations  are  shown  in  Table  l,  with  Pra  =  200  lb/ft2.  There 
are  several  reasons  for  this.  The  ARC  data  involved  a  propellant  for  which 
the  details  of  composition  were  not  readily  available,  and  the  questions 
underlying  the  validity  of  the  data  made  it  appear  that  only  an  approximate 


Table  1:  PROGRAM  INPUTS:  REFERENCE  INPUTS 


VEHICLE  PARAMETERS: 


FIRST  HALF  CONE  ANGLE 

THETA1  = 

6.00 

DEGREES 

SECOND  HALF  CONE  ANGLE 

THETA2  = 

.01 

DEGREES 

LENGTH  OF  VEHICLE 

XBODY  = 

3.75 

FEET 

LENGTH  OF  CONE 

YCONE  = 

2.75 

FEET 

THRUSTER  ANGLE 

THETAJ  = 

.00 

DEGREES 

FLIGHT  PARAMETERS: 

TEMPERATURE 

TINF 

=  390.00 

DEGREES  RANKINE 

ANGLE  OF  ATTACK 

ALPA 

=  0.00 

DEGREES 

BANK  ANGLE 

PHI 

=  0.00 

DEGREES 

VEHICLE  CENTER  OF  GRAVITY 

ICG 

=  1 .83 

FEET  (MEASURED  FROM  NOSE) 

RATIO  OF  SPECIFIC  HEATS 

GAMA 

=  1.40 

(DIMENSIONLESS) 

INJECTOR  PARAMETERS: 

CHAMBER  PRESSURE 

POJET  = 

50.00 

POUNDS  PER  SQUARE  INCH 
ABSOLUTE 

THROAT  AREA 

AO JET  = 

.03 

SQUARE  INCHES 

EXPANSION  RATIO 

ERATIO  = 

2.00 

(DIMENSIONLESS) 

INCLINATION  OF  JET 

ALFAJ  = 

-51 .00 

DEGREES 

INJECTOR  SITE  DISTANCE 

ZJET 

1  .60 

INCHES  (MEASURED  FROM 
VEHICLE  BASE) 

HALF  ANGLE  OF  NOZZLE  EXIT  CONE 

BETA 

15.00 

DEGREES 

PROPELLANT: 

FUEL 

NUMBER 

19 

ARCADENE  168(ARC 

TESTS) 

OXIDIZER 

NUMBER 

19 

ARCADENE  168(ARC 

TESTS) 

OXIDIZER-FUEL  RATIO 

OOVERF  = 

1.00 

(DIMENSIONLESS) 

PROGRAM  FLOW  PARAMETERS  : 

FROZEN/EQUILIBRIUM  CHEMISTRY 

ETAR 

_ 

0.00  (DIMENSIONLESS) 

NOSE  RADIUS  MULTIPLIER 

PCTR 

= 

.8000  (DIMENSIONLESS) 
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comparison  would  be  justified.  Because  of  the  essential  nature  of  the  injec- 
tant  annulus  In  the  calculation,  no  runs  without  injection  can  be  made,  and 
thus  only  a  qualitative  comparison  with  the  no-injection  data  can  be  made. 
Considering  these  limitations,  the  agreement  shown  in  Figure  12  is  believed 
to  be  reasonably  good. 


SECTION  FIVE 


PARAMETRIC  SURVEY 

Parameters  Studied 

A  large  number  of  inputs  are  required  for  the  calculations  described 
above.  Most  of  these  are  shown  in  Table  1,  which  also  shows  the  values  of 
these  inputs  that  are  taken  as  reference  conditions  for  the  parametric  survey. 
Table  2  shows  the  parameters  to  be  varied,  and  the  ranges  desired. 

From  Table  1  and  Figure  13,  the  vehicle  is  seen  to  be  a  biconic,  3.7b  ft. 
long  overall.  The  nose  is  a  6  deg.  half-angle  cone,  followed  by  an  almost- 
cylindrical  section  (0.01  deg.  half-angle),  1  ft.  long.  (The  program  cannot 
use  a  body  half-angle  of  zero.) 

The  atmosphere  is  assumed,  for  purposes  of  this  study,  to  have  a  constant 
temperature  of  390  R. 

Effect  of  Parameters 

Mach  Number  and  Altitude 

The  effect  of  flight  Mach  number  is  shown  in  Fiqure  14.  For  a  range  of 
altitudes,  the  continuous  variation  of  Pu/P  with  II  is  a  result  of  enforcing 
transition  between  plus  and  minus  solutions  as  described  above.  The  mixing 
half-angle  required  for  this  transition  was  o  =  1.0  deg.  This  is  an 
eminently  reasonable  value  —  an  incompressible  shear  layer  with  one  side 
bounded  by  stagnant  fluid  has  o  'v  6  deg.,  and  both  M  >  0  and  U2  >  0  (sub  2 
denoting  the  low-speed  side)  are  known  to  decrease  o. 


Table  2.  PARAMETERS  TO  BE  VARIED  (INCLUDING  DESIRED  RANGES) 

1.  Local  Mach  number  (2  to  5;  lower  limit  of  1.5 
if  possible) 

2.  Altitude  (sea  level  to  100K  ft) 

3.  Propellant  mass  flowrate 

4.  Missile  diameter  (base  diameter  from  2  to  16  in) 

5.  Propellant  type  (solids  with  and  without  A1 ; 
liquid  bipropellant;  monopropellant  liquid;  GH2) 

6.  Injector  location  on  missile  body 

7.  Injection  momentum 

8.  Number  of  injectors  (from  2  to  a  number  sufficient 
to  approximate  an  annulus) 


Table  2  shows  an  altitude  range  of  sea  level  to  100,000  ft  is  desired. 

In  terms  more  consonant  with  the  inputs  of  the  present  program,  this  corresponds 

2 

to  ambient  pressures  in  the  approximate  range  2000  <  P^  <■  20  lb/ft  .  Calcu¬ 
lations  cannot  be  made  for  a  single  vehicle  with  fixed  injection  conditions 
flying  over  this  range.  This  is  because  at  relatively  high  values  of 
("relative"  refers  to  the  jet  total  pressure  P  .)  the  injectant  flow  will 
separate  within  the  injectant  nozzle,  while  at  low  Pro  the  individual  jets 
produce  such  large  disturbances  that  the  jets  are  not  turned  downstream, 
even  approximately,  by  the  time  the  base  plane  is  readied. 

Thus  the  desired  altitude  range  can  only  ue  covered  by  varying  Pqj  with 
Pa,  so  as  to  stay  within  ttie  range  for  which  calculations  can  be  made.  Because 
of  the  importance  of  altitude  as  a  parameter,  these  calculations  are  made 
over  the  entire  Mach  number  range.  Results  are  shown  in  Figure  14. 

Several  aspects  of  these  results  are  of  interest.  First,  it  appears  that 

the  case  chosen  to  show  the  effect  of  Mach  number  in  Figure  12,  the  choice 

2 

being  made  because  Pm  =  200  lb/ft  corresponds  to  an  altitude  intermediate 
in  the  desired  range,  was  somewhat  unique  in  that  performance  at  all  other 
pressures  is  increasing  with  Mach  number  at  =  2.  It  is  not  clear  whetner 
this  trend  is  a  real  one,  or  an  artifact  of  the  approximations  used. 

Second,  performance  is  highest  at  the  highest  altitude.  This  does  not 

include  combustion,  and  since  P  .  is  scaled  with  P  we  might  expect 

oj  00 

identical  performance  to  be  obtained  here  (and  in  fact  over  the  whole  range). 
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The  difference  Is  believed  to  arise  from  the  use  of  fully  equilibrium  chemistry 
to  obtain  jet  conditions  in  the  individual-jet  calculation. 

Injectant  Mass  Flow  Rate 

This  important  parameter  can  be  varied  in  more  than  one  way.  Two  of  the 
simplest  were  chosen  for  the  present  survey:  variation  of  jet  throat  area 
Aojet  and  Pressure  p0jet*  Resu^ts  varying  these  inputs  are 

shown  in  Figure  lb,  for  =  2  and  o  *  l.U  deg,  and  for  the  reference 
conditions  of  Table  1. 

Several  interesting  things  can  be  seen  from  this  Figure.  Variation  of 
m.  by  the  two  methods  has  opposite  effects.  This  is  consistent  with  existing 

J 

data,  which  indicate  that  low-velocity  injection  is  most  effective  for  a 
14  29 

unit  of  mass  flow,  *  and  which  may  be  interpreted  as  meaning  that  the 
best  performance  is  reached  by  disturbing  the  external  stream  least. 

•  • 

The  behavior  of  PR/P  as  m.  0  is  of  interest.  The  data  range  for  m.  =  0 

o  00  j  j 

is  shown  on  the  ordinate  in  Figure  15.  As  AQjet  ->•  0  the  computed  P^/P^ 
is  clearly  approaching  a  value  that  will  at  least  approximately  agree  with 
experiment;  just  as  clearly,  when  PQjet  -+0.  the  correct  limiting  value  will 
be  reached  only  by  what  seems  to  be  an  improbable  change  in  shape  of  the  curve. 
This  behavior  is  so  far  unexplained  and  needs  an  explanation. 

Finally,  Figure  15  also  shows  a  large  scatter.  All  points  shown  are  actual 
calculations.  This  sort  of  scatter  did  not  occur  in  runs  over  the  Mach  number 
or  altitude  ranges  previously  discussed,  and  to  find  it  here  is  an  unpleasant 
surprise.  Its  cause  is  not  clear,  numerical  solutions  can  be  expected  to 
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show  scatter  of  the  order  of  the  convergence  criteria;  but  here  G.T*  is  used, 
and  we  observe  scatter  of  order  10".  Of  course  a  number  of  such  criteria 
must  be  satisfied,  but  even  neglecting  their  almost  certain  independence  and 
summing  their  effects  can  account  for  less  than  an  order  of  magnitude  increase 
in  scatter,  leaving  another  order  of  magnitude  to  he  accounted  for.  (The 
scatter  caused  by  a  number  of  truly  independent  sources  is  found  as  the 
square  root  of  the  sum  of  the  squares  of  the  individual  values,  which  is  less 
than  their  simple  sum,  as  the  hypotenuse  of  a  right  triangle  is  shorter 
than  the  sum  of  the  other  two  sides.) 

(Jne  possibility,  which  has  not  as  yet  been  made  definite,  is  that  part  of 

the  problem  lies  in  the  chemistry  subroutines.  In  the  original  development  of 

dBLIP  it  was  found  that  the  chemistry  was  very  slow  to  equilibrate,  in  the 

sense  of  reaching  a  constant  value  of  Ti+^,  even  for  frozen  flow;  and  that 

(R+ll  (kl 

even  after  a  criterion  of  the  form  e  =  1  -  T.+j  VTi+i  '  was  satisfied, 

slow  changes  in  T^+,  with  increasing  k  (the  loop  number)  were  observed.  Even 

in  the  present  formulation,  after  satisfying  <  0.001,  an  arbitrarily  chosen 

100  loops  through  the  conservation  equations  and  thermo  chemistry  are  made 

to  ensure  equilibration.  It  may  be  that  even  with  this  provision,  a  fully 

consistent  value  of  T..,  has  not  been  found. 

i  +  l 

Missile  Diameter 

The  reference  vehicle  (Table  1  and  Figure  13)  is  a  biconic  that 
approximates  a  cone-cylinder;  the  nose  is  a  b  deg.  half-angle  cone,  Z.  7s  ft. 
(33  in.)  long,  and  the  "cylindrical"  section  is  actually  a  0.01  deg.  half¬ 
angle  cone,  one  foot  long.  (The  "cylindrical"  section  thus  diverges  by 
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U.GU2U9  in.  in  radius.)  The  reference  base  diameter  is  b.93b  in.  To  vary 
the  base  diameter  over  the  desired  range  we  vary  the  cone  length  over  the 
range  (in  feet)  0.833  <  XCUHL  <  b.67. 

If  all  lengths  in  the  problem  are  scaled  to  keep  geometrical  similarity, 
the  same  base  pressure  would  result.  Tnis  is  also  approximately  true  ?f 
the  “cylindrical"  section  length  is  held  fixed  as  other  lengths  chanqe,  since 
although  local  conditions  are  affected  by  the  proximity  of  the  conical  section, 
tnis  effect  is  not  large.  We  therefore  take  the  more  difficult  route  of 
holding  all  other  variaoles  constant  while  XCGht  is  varied.  Wnen  this  is  done. 
Figure  lo  results.  It  shows  base  pressure  increasing  with  decreasing  Lase 
diameter.  This  can  be  viewed  as  equivalent  to  an  increase  in  the  nondimensional 
injection  flowrate  nij/p^  v^  Ag;  as  shown  in  Figure  lu,  a  base  diameter  of  2  in. 
yields,  with  tne  reference  conditions  of  Table  1,  a  nondimensional  injection 
rate  of  0.207.  This  method  of  changing  flowrate  may  be  considered  equivalent 
to  changing  A  -et,  since  tne  jet  velocity  is  not  changed  in  either  case 
(see  the  discussion  above).  If  so,  we  can  combine  the  data  shown  in  Figures 
lb  and  lb  to  yield  tne  result  shown  in  Figure  17. 

iio  line  is  drawn  through  the  different  regions  in  Figure  17,  as  tney  were 
obtained  by  different  methods  and  it  is  not  certain  that  comparison  is  valid, 
nevertheless ,  if  one  were  known,  its  snape  would  be  quite  reasonable.  Its 
slope  would  be  proportional  to  specific  impulse,  and  would  clearly  be  greater 
at  low  than  at  high  injection  rate.  This  is  consistent  with  all  experience 
in  tiiis  and  similar  areas.  It  would  apparently  require  very  substantial 
injection  to  raise  the  base  pressure  auove  ambient,  and  this  too  is  consistent 
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Figure  17.  Effects  of  Varying  Nondimensional  Injection  Rate  by  Different  Means 


with  expectations. 


Propellant  Type 

Runs  were  made  for  the  reference  propellant,  an  aluminized  solid  propellant, 
as  well  as  for  a  representative  gas  generator  (non-aluminized)  solid  and  a 
liquid  bipropellant.  Properties  of  each  injectant  are  shown  in  Tables  3, 

4  and  5  respectively;  compositions  of  the  exhaust  gases  in  the  combustion 
chambers  are  shown  respectively  in  Tables  6,  7  and  8. 

Results  are  shown  in  Table  9.  The  aluminized  solid  is  significantly  better 
than  the  non-aluminized  solid,  and  it  in  turn  is  better  than  the  liquid 
bipropellant.  The  latter,  however,  is  run  at  the  approximately  stoichiometric 
0/F  ratio  of  one;  a  lower  value  might  well  aid  performance. 

Runs  could  not  be  made  with  heated  gaseous  hydrogen  (GH2)  nor  for  the 
liquid  monopropellant,  although  it  appears  the  problems  encountered  are  rela¬ 
tively  trivial.  The  GH^  case  will  probably  run  if  the  correct  shock  scale  can 
be  found;  the  work  done  finding  this  value  for  more  conventional  propellants 

does  not  apply  directly  to  hydrogen,  as  substantial  molecular  weight  effects 
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are  likely  and  in  fact  have  been  observed.  The  information  required  to  repeat 
this  work  is  believed  to  be  available.  The  liquid  monopropellant  to  be  used 
is  hydrazine,  whose  properties  are  available  in  the  program  as  a  liquid 
injectant  only.  Inclusion  of  its  gaseous  products  of  decomposition,  H2  and  N2, 
will  very  likely  allow  runs  to  be  made. 


Table  3.  REFERENCE  PROPELLANT 


JET  CONDITIONS: 


MACH  NUMBER 

MJET 

1.9781 

(DIMENSIONLESS) 

EXIT  TEMPERATURE 

TJET 

2313.68 

DEGREES  RANKINE 

EXIT  VELOCITY 

VJET 

5587.27 

FEET/SECOND 

EXIT  DENSITY 

ROJET  = 

.0001396 

SLUGS/CUBIC  FOOT 

MOLECULAR  WEIGHT 

MOLWTJ  = 

31.047 

(DIMENSIONLESS) 

RATIO  OF  SPECIFIC  HEATS 

GAMAJ  = 

1.1963 

(DIMENSIONLESS) 

JET  THRUST 

FR 

11.48 

POUNDS 

JET  MASS  FLOW  RATE 

MDOTJ  = 

.06 

POUNDS/SECOND 

JET  THERMOCHEMISTRY: 

TOTAL  ENTHALPY 

PROPELLANT 

-362.620 

CAL/ GRAM 

PRODUCTS 

-1434.603 

CAL/GRAM 

HEAT  RELEASE 

1091.983 

CAL/ GRAM 

EXIT  ENTHALPY 

-709.037 

CAL/GRAM 

SPECIFIC  HEAT 

.39012 

CAL/GRAM/DEG 
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Table  4.  NON-ALUMINIZED  SOLID  PROPELLANT 


JET  CONDITIONS: 


MACH  NUMBER 

MJET 

1.9991 

(DIMENSIONLESS) 

EXIT  TEMPERATURE 

TJET 

971.25 

DEGREES  RANKINE 

EXIT  VELOCITY 

VJET 

4916.93 

FEET/SECOND 

EXIT  DENSITY 

ROJET  = 

.0001839 

SLUGS/CUBIC  FOOT 

MOLECULAR  WEIGHT 

MOLWTJ  = 

19.057 

(DIMENSIONLESS) 

RATIO  OF  SPECIFIC  HEATS 

GAMAJ  = 

1.3264 

(DIMENSIONLESS) 

JET  THRUST 

FR 

1.95 

POUNDS 

JET  MASS  FLOW  RATE 

MDOTJ  = 

.01111 

POUNDS/SECOND 

JET  THERMOCHEMISTRY: 

TOTAL  ENTHALPY 

PROPELLANT 

-424.480 

CAL /GRAM 

PRODUCTS 

-961.359 

CAL/GRAM 

HEAT  RELEASE 

536.879 

CAL/GRAM 

EXIT  ENTHALPY 

-692.760 

CAL/GRAM 

SPECIFIC  HEAT 

.42377 

CAL/GRAM/DEG 
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Table  5.  LIQUID  BIPROPELLANT:  NITROGEN  TETROXIDE/ 
MONOMETHYL  HYDRAZINE  (O/F  =  1.0) 


JET  CONDITIONS: 


MACH  NUMBER 
EXIT  TEMPERATURE 
EXIT  VELOCITY 
EXIT  DENSITY 
MOLECULAR  WEIGHT 
RATIO  OF  SPECIFIC  HEATS 
JET  THRUST 
JET  MASS  FLOW  RATE 


MJET 

= 

2.0952 

TJET 

= 

1464.64 

VJET 

= 

6677.35 

ROJET 

= 

.0001004 

MOLWTJ 

= 

16.752 

GAMAJ 

= 

1.2981 

FR 

= 

1.94 

MDOTJ 

= 

.00823 

JET  THERMOCHEMISTRY 


TOTAL  ENTHALPY 
PROPELLANT 
PRODUCTS 
HEAT  RELEASE 
EXIT  ENTHALPY 
SPECIFIC  HEAT 


114.500 

-931.425 

1045.925 

-380.276 

.51649 
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(DIMENSIONLESS) 

DEGREES  RANKINE 

FEET/SECOND 

SLUGS/CUBIC  FOOT 

(DIMENSIONLESS) 

(DIMENSIONLESS) 

POUNDS 

POUNDS/SECOND 


CAL/GRAM 
CAL/GRAM 
CAL/GRAM 
CAL/GRAM 
CAL/GRAM/DEG  K 


Table  6.  REFERENCE  PROPELLANT:  CHEMICAL  COMPOSITION  OF  PRODUCTS 


SPECIES 

MOLE  PERCENT 

WEIGHT  PERCE 

H2 

48.4507 

3.1461 

CO 

24.8259 

22.3980 

ALCL 

13.2544 

26.6549 

AL203* 

9.6359 

31.6444 

N2 

7.1218 

6.4264 

AL 

2.6498 

2.3027 

AL20 

1.8733 

4.2211 

H 

1.0867 

.0353 

AL203$ 

. 549926E+00 

. 1 80597E+01 

ALCL2 

. 279025E+00 

.879782E+00 

HCL 

.250906F+00 

. 294688E+00 

ALH 

. 160925E+00 

. 145068E+00 

CNH 

. 300251 E-01 

. 261 372E-01 

H20 

. 769029E-02 

. 446249E-02 

CL 

.405249E-02 

.462807E-02 

ALOCL 

. 302404E-02 

. 763983E-02 

C02 

.684904E-03 

.970883E-03 

ALO 

.406398E-03 

. 562592E-03 

C2H2 

.  20101 7E-03 

. 168584E-03 

ALCL3 

.  1 66592E-03 

. 71 5529E-03 

CHO 

.113141E-03 

.105750E-03 

CH4 

.  600399E-04 

. 310243E-04 

CH3 

.  485094E-04 

. 23491 2E-04 

NH3 

.418640E-04 

. 229658E-04 

HO 

. 290862E-04 

. 1 59337E-04 

AL202 

. 218865E-04 

. 605966E-04 

CN 

. 938452E-05 

. 786463E-05 

NH 

. 673214E-05 

. 325599E-05 

NH2 

. 546039E-05 

. 281 81 9E- 05 

ALN 

. 369861 E-05 

. 488282E-05 

CH20 

. 324710E-05 

. 31 4039E-05 

N 

. 1 78230E-05 

.804142E-06 

ALHO 

. 779190E-06 

. 1 10396E-05 

NO 

. 585203E-06 

. 56561 3E-06 

0 

. 425825E-06 

. 219446E-06 

ALC 

. 201 274E-06 

. 262771 E-06 

CH2 

. 1 34244E-06 

.606505E-07 

C2N2 

. 109581 E-06 

. 1 83667E-06 

C 

. 710821 E-07 

. 274989E-07 

C2H4 

.679077E-07 

. 61 3606E-07 

CH 

. 678209E-07 

. 284392E-07 

CL2 

. 164527E-07 

. 375790E-07 

Table  7.  NON-ALUMINIZED  SOLID  PROPELLANT 


CHEMICAL  COMPOSITION  OF  PRODUCTS 


SPECIES 

MOLE  PERCENT 

WEIGHT  PERC 

H2 

39.6783 

4.1974 

CO 

33.1352 

48.7028 

HCL 

11.3475 

21.7126 

N2 

6.1972 

9.1105 

C02 

5.0540 

11 .6715 

H20 

3.6153 

3.4178 

CH4 

. 970899E+00 

.81 71 58E+00 

C$ 

. 584828E+00 

. 368590E+00 

NH3 

. 167865E-02 

. 1 50024E-02 

CNH 

. 859082E-04 

. 1 21 827E-03 

CH20 

. 757954E-05 

. 1 1 9424E-04 

CH3CL 

. 620245E-05 

. 1 64332E-04 

C2H4 

.3494868-05 

. 514470E-05 

CH3 

. 1 59218E-06 

. 12561 2E-06 

H 

.  103375E-06 

. 546782E-08 

CL 

. 28231 4E-07 

. 525255E-07 

CHO 

. 24811  IE-07 

. 377802E-07 

C2H2 

. 232478E-07 

. 31 7633E-07 
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Table  8.  LIQUID  BIPROPELLANT:  NITROGEN  TETROXIDE/ 
MONOMETHYL  HYDRAZINE  (O/F  =  1.0) 


CHEMICAL  COMPOSITION  OF  PRODUCTS 


SPECIES 

MOLE  PERCENT 

WEIGHT  PERCE 

H2 

38.8847 

4.6795 

N2 

27.2821 

45.6258 

H20 

15.6523 

16.8331 

CO 

15.5993 

26.0832 

C02 

2.5800 

6.7781 

H 

.1 16831 E-02 

. 702982E-04 

NH3 

. 290985E-03 

. 295843E-03 

CNH 

. 1 31 782E-04 

. 21 2575E-04 

HO 

. 1 21 523E-04 

. 1 23378E-04 

CH4 

. 66291 6E-05 

.634849E-05 

CH20 

. 1 91991 E-05 

. 344127E-05 

CHO 

. 162026E-05 

.280669E-05 

NO 

. 106044E-06 

. 189955E-06 

NH2 

.61 937  3E-07 

. 592446E-07 

i 


I 

1 


i 


Table  9.  EFFECT  OF  PROPELLANT  TYPE 


M  =2.5 

'JO 

a  =1.0 

REFERENCE  CONDITIONS 

Propellant  Type 
Reference  Solid 

Non-Aluminized  Sol  id 

Liquid  Bipropellant  (NTO/MMH,  0/F  =1.0) 


0.67 

0.60 

0.51 
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Injector  Location 

The  streamwise  location  of  injection  nozzles  is  a  likely  field  for 
optimization.  The  present  program  cannot  treat  nozzles  nearer  the  base 
than  1.33  h,  where  it  is  ttie  jet  "penetration;"  nor  can  it  treat  nozzles  so 
far  upstream  that  the  individual-jet  shock  is  too  close  to  the  Mach  angle. 

Tne  former  is  a  fundamental  limitation  (although  the  factor  1.33  might  oe 
decreased  toward,  but  never  below,  one);  tne  latter  is  not,  and  might 
possibly  be  removed.  Within  tiiese  limits,  calculations  have  been  made  varying 
the  distance  on  the  vehicle  surface  from  the  nozzle  centerline  to  tne  ease 
plane,  z0jet-  Results  are  shown  in  Figure  18.  They  are  in  accord  with 
expectations. 

Injector  Momentum 

Performance  is  calculated  for  varying  values  of  jet  velocity  at  fixed  jet 

flow  rate.  The  simplest  way  of  doing  this  is  oy  changinq  the  injection  nozzle 

area  ratio.  Results  of  doing  so  are  shown  in  Figure  19.  The  decrease 

of  base  pressure  with  increasing  jet  velocity  is  at  least  consistent  with 

29 

results  of  Schadow  and  Chieze,  although  direct  comparison  is  not  possible 
between  axisymmetric  calculations  and  planar  flow  experiments. 

dumber  of  Injectors 

A  correction  has  been  included  in  the  program,  almost  since  its  earlier,; 

existence,  for  the  failure  of  the  individual  jet  flows  to  completely  surround 
14 

the  base  cavity.  It  is  based  on  limited  data  and  on  the  assumption  t.iat 
tne  correction  is  of  the  form 

Js STL  _ gap  width _ 

Pp  H  base  circumference 

uncorr 
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0  12  3 

Zojet(IN) 

Figure  18.  Effect  of  Injector  Location 


63 


NOZZLE  AREA  RATIO  ERATIO 

i _ i  ■ 

4,991  5,587  6,148 

JET  VELOCITY  VJET  (FT/SEC) 

Figure  19.  Effect  of  Injectant  Momentum  Flux  at  Fixed  Mass  Flowrate 
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where  the  gap  width  is  the  total  length  of  the  base  perimeter  that  is 

»  not  covered  by  the  inner  flows  calculated  for  the  individual  jets,  and 

constants  multiplying  the  exponential  function  and  its  argument  are  found 

14 

from  the  existing  data  points,  of  winch  there  are  two. 

In  addition  to  this  correction,  there  are  ottier  effects  of  the  number 
of  nozzles  A  which  appear  in  the  calculation.  As  N  is  decreased  the  scale 
of  the  individual  jets  increases  with  respect  to  the  vehicle;  the  base  plane 
conditions  (the  annulus  entrance  conditions)  consequently  change,  both 
because  the  end  of  the  turning  region  is  shifted  closer  to  the  base  plane, 
and  because  the  body  transverse  curvature  may  become  important  (jet 
wraparound).  These  effects  are  shown  in  Figure  20.  The  desired  range  of 
N  could  not  be  covered  because  for  N  <  8  the  individual-jet  turning  regions 

I 

come  too  near  the  base.  In  the  absence  of  data,  the  predictions  shown  in 
Figure  20  are  in  accord  with  the  expectation  that  as  11  increases  performance 
should  also  increase. 

Other  Parameters 

In  developing  this  parametric  survey,  emphasis  was  of  course  given  to 
those  parameters  that  were  explicitly  requested,  liany  other  parameters  are  of 
interest  and  can  be  accounted  for  in  the  present  program.  Here  a  few  of  these 
that  seem  important  are  discussed. 

.  Combustion  Efficiency;  This  variable,  denoted  n(,  or  ETAR,  represents  the 

fraction  of  the  annulus  flow  (or  in  the  indi vidual-jet  calculations,  of  the 
inner  flow)  that  is  treated  using  chemical  equilibrium.  Based  on  the  only 

i 
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Figure  20.  Effect  of  Number  of  Nozzles  ,  Fixed  I  njectant  Flow  Rate 


14 

conclusive  data  on  the  subject,  the  value  =  U  was  used  for  all  runs 
shown  in  this  report.  Runs  for  >  0  can  be  made,  and  will  be,  in  order 
to  better  understand  the  operation  of  the  program. 

boattailinq:  It  has  been  observed  that  good  performance  can  be  achieved  by 
locating  the  jets  on  a  boattail.  Such  a  geometry  can  probaoly  be  approximated 
in  the  present  program,  though  some  changes,  probably  small,  may  be  needed  to 
handle  a  negative  0^.  Tne  limitation  to  a  biconic  vehicle  means  that  the  body 
upstream  of  the  boattail  cannot  oe  a  true  cylinder,  and  that  even  if  6^  were 
made  small,  the  body  length  would  be  extreme.  This  is  not  a  serious  problem, 
however,  as  tiie  length  has  only  a  small  effect  on  the  inviscid  flow. 

nozzle  Inclination:  This  variable  (“jet)  determines  the  distribution  of  overall 
thrust  between  the  axial  component  of  nozzle  thrust  and  the  rise  in  base  pressure. 
As  such  it  is  a  natural  candidate  for  optimization  studies.  It  can  easily  be 
changed  in  the  present  program,  although  a  better  understanding  of  the 
einpi ri cal -shock  method  may  be  needed  to  run  at  the  lowest  Mach  numbers. 
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SECTION  SIX 


CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

The  results  presented  above  to  a  large  extent  cover  the  desired  ranges  of 
variables,  but  they  do  not  do  so  completely.  The  main  failure  is  in  treating 
some  injectants,  although  it  is  believed  that  the  problems  in  doing  so  are  not 
prohibitive.  Program  achievements  include  coverage  of  a  range  of  Mach  numbers 
that  was  doubtful  at  the  beginning  of  the  work,  of  two  orders  of  magnitude  in 
ambient  pressure,  and  of  substantial  ranges  of  other  parameters  of  interest. 
Some  uncertainty  is  present  in  the  results,  particularly  with  respect 
to  the  scatter  that  sometimes  results  (e.g.  Figure  lb),  as  well  as  in  the 
method  chosen  to  eliminate  multiplicity  of  solutions.  Trends  and  values 
calculated,  however,  generally  match  data  and  expectations  reasonably  well. 

Recommendations  for  Future  Work 

Although  many  of  its  results  compare  reasonably  well  with  data  and 
expectations,  the  model  and  associated  code  described  herein  fall  short  of 
what  is  desired  in  several  ways.  Symptoms  include  the  inability  of  run 
with  certain  propellants  (a  difficulty  which  is  probably  not  fundamental  a 
probauly  can  be  cured  with  relative  ease);  and  the  more  fundamental  problems  - 
absent  or  multiple  solutions,  scatter,  and  the  few  questionable  trends  (e.g. 
the  behavior  of  Pu/P  as  P  .  -*-U;  the  change  in  slope  of  Pn/P  vs.  M  at  M 

D  oo  OJ  D  00  m 

There  are  also  areas  where  the  physical  representation  is  inadequate.  The 
most  important  of  these  is  the  failure  to  include  an  inner  shear  layer  and 
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the  resulting  viscous-inviscid  interaction  downstream  of  reattachment. 

Failure  to  explicitly  satisfy  the  radial  momentum  equation  is  in  principle 
equally  important,  although  as  a  practical  matter  the  error  involved  is  in 
all  cases  observed  to  be  relatively  small  (<  17%)  and  to  probably  have  only 
a  negligible  effect  on  trends.  Other  physical  approximations  that  may  be 
of  lesser  importance  include  the  assumptions  of  zero  radial  pressure  gradient, 
of  single-ray  turning  of  the  injectant  annulus,  of  neglect  of  one  family 
of  characteristics  in  the  external  flow,  of  instantaneous  transition  from 
individual  jets  to  an  annulus  at  the  base  plane,  and  many  others. 

Of  these  shortcomings,  two  are  believed  most  in  need  of  additional  work. 
These  are 

(1)  A  deeper  understanding  of  the  problem  we  see  as 
multiple  solutions. 

(2)  Inclusion  of  the  inner  shear  layer  and  viscous-inviscid 
interaction. 

These  may  be  interrelated,  a  possibility  which  makes  tfie  second  of  higher 
priority.  As  shown  in  Figure  1 U  and  the  accompanying  discussion,  multiple 
solutions  (wriich  includes  its  alternate  possibility,  no  solution)  arise  uecause 
pressures  at  both  inner  and  outer  boundaries  of  the  injectant  annulus  are 
decreasing  functions  of  t)  ,  the  wake  cavity  half-angle.  Physically  it  is 
likely  that  the  inner  boundary  pressure  will  begin  to  rise  with  increasing 
6W  at  large  enough  ow.  This  is  because  at  very  large  0  (say  near  Du  deg.) 
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a  large  fraction  (perhaps  nearly  half)  of  the  annulus  flow  would  be  forced 
to  recirculate;  and  a  relatively  nigh  pressure  at  the  cavity  boundary  would 
ue  required  to  turn  inward  the  radially  outward  momentum  of  this  flow  at  and 
near  the  base  plane. 

If  the  inner  boundary  pressure  rises  with  increasing  8W  this 
guarantees  a  solution,  subject  only  to  tne  proviso  that,  at  some  smaller  ■ 

f’inner  <  Pouter*  T,lis  in  turn  cilan9es  t,ie  l'"ole  method  of  obtaining  solutions, 

and  in  addition  improves  solution  stability.  (The  latter  may  arise  from  the 

small-angle  crossings  shown  in  Figure  lb,  and  may  in  fact  oe  related  to  the 

scatter  in  some  of  the  program  predictions  that  is  referred  to  above. ) 

Tne  interrelation  comes  from  the  fact  that  item  (2)  above,  tne  inclusion  of 

the  inner  shear  layer  and  the  interaction  of  its  viscous  wake  with  the 

approximately  inviscid  annulus  and  external  flows,  is  required  in  order  to 

see  tuis  increase  of  P.  „  witn  9.  that  is  expected  on  physical  grounds. 

inner  w 

Although  the  work  described  above  should  probably  be  given  priority,  other 
work  is  also  recommended,  though  probably  at  a  lower  level  of  importance. 
Several  other  assumptions  have  been  mentioned  above,  and  they  are  not  a 
complete  list.  The  inaccuracies  resulting  are  probably  smaller  but  not 
negligible.  It  is  recommended  that  estimates  ue  made  of  the  magnitude  o'  t 
inaccuracy  generated  by  at  least  the  apparently  most  important  of  these 
assumptions  (including  those  judged  more  important,  as  above).  When  this  i  » 
done,  an  informed  assessment  both  of  the  expected  program  accuracy  and  the 
best  subject  for  further  work  can  ue  made. 
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NOMENCLATURE 


Area 

Vehicle  base  area 
Injectant  nozzle  throat  area 

Jet  scale  length  ("penetration  height");  altitude 

Local  Mach  number 
Flight  Mach  number 

Mass  flowrate 

Injectant  mass  flowrate 

(If  no  units  are  specified,  units  are 

Number  of  nozzles 

Pressure 

Radius 

Velocity 

Length  of  vehicle 
Length  of  vehicle  nose  cone 


Angle  of  equivalent  single  ray 

Inclination  of  jet  with  respect  to  normal  to  centerline; 
positive  upstream 

Ratio  of  specific  heats 

Annulus  divergence  angle 

Mass  entrained  downstream  of  expansion  ray 

Mass  entrained  upstream  of  expansion  ray 

Reaction  efficiency  (fraction  of  flow  at  equilibrium; 
remainder  frozen) 

Turning  angle 

First  cone  half-angle  of  biconic  vehicle 

Second  cone  half-angle  of  biconic  vehicle 

Second  cone  half-angle  of  biconic  vehicle  (=62) 

Cavity  half-angle 

Mach  angle 

Prandtl -Meyer  angle 

Density 

Mixing  half-angle 


SUBSCRIPTS 


1  Annulus  leading  edge 

i+1  Annulus  trailing  edge 

in  Annulus  inner  boundary 

out  u  Annulus  outer  boundary,  upstream  of  expansion  ray 

out  d  Annulus  outer  boundary,  downstream  of  expansion  ray 

1  External  flow  at  annulus  leading  edge 

2  External  flow  immediately  after  expansion  ray 

0  Stagnation  values 


, 
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